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rG . Abstract 

I , Discrete-time dynamical systems can be derived from group actions. 

rG ' In the present work possibility of application of this method to systems 

jrt I of ordinary differential equations is studied. Invertible group actions are 

considered as possible candidates for stroboscopic maps of ordinary differ- 
ential equations. It is shown that flow of the Bloch equation is a unique 
suspension of an invertible map on the SU (2) group. 
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Discrete-time dynamical systems can be formulated in terms of group actions to 
exploit the group structure and get a better understanding of the correspond- 

^T '_ ing dynamics. This approach was used to study non-invertible discrete-time 

dynamics on E (2) [1] and SU (2) [H [H [H [4]. The Shimizu-Lcutbecher map 
[51 [S], a tool to study group structure, was solved for an arbitrary Lie group 
G [I]. Since the Shimizu-Leutbecher sequence generates the logistic map for 
Q = SU (2) [U 13] a general solution of the logistic map was thus obtained. For 

^ . Q = E{2) this approach led to a better understanding of the Harter-Heighway 

$—( ' fractal curve [1] . A two-dimensional generalization of the logistic map was also 

introduced as a map on SU (2) and investigated [4|. 

On the other hand, structure of Kleinian groups is naturally studied in the 
setting of discrete-time dynamical systems, revealing in this way connections 
with fractals [3 [HI [S] • There were several attempts to study dynamics on Lie 
groups within continuous-time rather than discrete-time approach. Continuous- 
time dynamics on SU (2) and SU (2) x SU (2) was defined in the setting of 
ordinary differential equations (ODEs) [THl [H] while continuous iteration of 
maps was defined to find a correspondence between mappings and continuous- 
time evolution [TH [I^ • 

The question now arises whether a general connection between discrete-time 
group actions and continuous-time dynamical systems can be established. 



Wc investigate a possibility of relating group actions with stroboscopic maps 
of ODEs. Let us consider a continuous-time dynamical system given by a set of 
ODEs 

^=f(x(t),t), (1) 

where x = [xi, X2, ■ ■ ■ , a;„] , f = [/i, /2, ... , /„]. Let x {t) be a solution to 
Eq.([T]). Then the invertible map 

St. X (t) ^— . X (t + T) , (2) 

is a stroboscopic map of strobe time T. Stroboscopic maps with several strobe 
times are standard tools to solve differential equations to mention only the 
Runge-Kutta methods. On the other hand, a discrete map can be also used 
to generate a continuous flow. Such a flow, non unique in general, is called a 
suspension of the map [141 HI] . 

In the present work we shall consider a class of invertible maps on a group 
and then try to find a flow, preferably unique, such that its trajectories connect 
continuously the iterates of the map. 

The paper is organized as follows. In the next Section maps on a Lie group 
Q are defined and method to deduce evolution of group parameters is described. 
Principal results are described in Sections 3 and 4. In Section 3 a class of 
invertible maps is considered. These maps are solved in Subsection 3.1. In 
Subsection 3.2 a simplified map is defined, solved and parameterized on the 
SU (2) group. In Subsection 3.3 symmetry and restrictions of the dynamics are 
determined. Results of the Subsection 3.2 are used in Section 4 to derive ODE 
for which the simplified map is a stroboscopic map. The map samples the fiow 
of the ODE exactly and arbitrarily densely and it follows that the construction 
is unique. This ODE is the Bloch equation [l6j, also known as the Landau- 
Lifshitz equation. These results generalize an earlier findings |17]. In Section 
5 computations for the simplified map on the SU (2) group are presented to 
elucidate dynamics of the Bloch equation. In the last Section the obtained 
results are summarized and perspectives of further research are outlined. 

2 Group dynamical systems 

Let us consider a dynamical system defined by the following map 

Gn+1 = fiGN, . . .) , (3) 

where Gn G G- Let Q he a, simple Lie group and g its Lie algebra. Then Gn G G 
can be written in exponential form 

Gn = exp (Xn) , (4) 

where X^ G |18| . Any infinitesimal operator X^ of a n-parameter Lie group 
Q is a. linear combination of n generators I^ 

n 

X^ = ^/^4, I^es, (5) 

fc=i 



where real parameters c^, ..., c^ are local coordinates of the Lie group element 

Substituting Eqs.Q, ^ into Eq.Q and using completeness of the basis 
consisting of the generators /'' and the unit matrix 1 we get a discrete-time 
dynamical system in parameter space [TJ 119] 

cjv+i =^-' (^iv> •■•. c^) > J ^'^, ■■■, n, (6) 

where F^ are continuous functions pp] . 

3 Discrete-time dynamics on a group 

Let us consider an invertible discrete-time dynamical system on a Lie group Q 

Rn+1 = QnRnQn-iRn-iQn^iR^j Qpj , N — 1,2,..., (7) 

i.e. (fi = QnRnQn-iRn-iQn'-iRn^Qn^ in Eq.(I3|), where we assume knowl- 
edge of all group elements Qm needed for the computations. Let us note 
here that apparently simpler non- invertible Shimizu-Leutbecher map Rn+i = 
RnQRJ;^ has a complicated solution [1]. It turns out, however, that a relatively 
simple solution to ([7]) can be constructed upon introducing new variables Sn 
defined below. 

3.1 Exact solution 

We note that Eq.([7]) can be reduced to two simpler equations. Indeed, intro- 
ducing new quantity 

Sn — RnQn-iRn-iQn-2, (8) 

we can rewrite Eq.([7|) as 

Rn+i ~ Sn+iRn-iSj^^i, N ~ 1,2, . . . , (9a) 

Sn+1 - QnSnQn'-2^ N=l,2,... . (9b) 

To run dynamics defined by ([7]) or, alternatively, by ^, we have to impose 
initial conditions for Eq.(I7]), Rq, Ri, from which initial condition for Eq. (|9bD . 
^2 = R2Q1R1Q0, can be also computed. Equations ^ are easily solved 



R2K — S2KS2K-I ■ ■ ■ S2R0S2 ■ ■ ■ S2K~lS2K^ (10a) 
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Sn — Qn-1 ■ ■ ■ QiSiQ_i . . . Qn-3j (10b) 



where K — 1,2,... , N = 2, 3, . . . and similar equations can be written for 

R2K+1- 



3.2 Discrete-time dynamics on the sphere 

In the case Q = SU{2) the Hamilton's parameterization can be used. Let unit 
vectors rjv = {rjf, r%, r^),q7v = {qn, q%, q%) and angles x at, aw correspond 
to rotation axes and rotation angles, respectively. Then the rotation matrices 
Rn, Qn are defined as 

Rn = exp(i^aTw), |rAr| = l, (11) 

Qn = exp (i^CT • qw) , IqwHl, (12) 

where i^ ~ —1, and a — [a^, a^, a^] is the pseudo vector with the Pauli 
matrices as components [H] . With parameterization pT]) , ([TH ) Eq. ([7]) induces 
dynamics of vectors r^ on unit sphere. We shall consider a special case Qn = Q 

Rn+1 = QRnQRn-iQ Rn Q ■ (1^) 

The solution of (fT5|) is immediately obtained from PH)) : 

R2K - Q'^^'P^'R^P-'^Q-^'', (14a) 

R2K+1 = Q2^P^i?iP^Q-2^, (14b) 

where P i Q-^SiQ-^ = Q-^RiQR^, R2K+1 = Q-^R2K+iQ, X = 0, 1, 2, ... . 
It follows that equations (|14a|) . (J14bp generate analogous dynamics. Matrices 
Q, P are parameterized as 

Q = exp (j|cr • q) , |q| = 1, (15a) 

P^expfzfa-pV |p| = l, (15b) 



where a, (3 are the corresponding angles of rotations while p = [p"'^,p^,p"^] and 
q = [9"'^7 9^i9^] are unit vectors. We still have to impose initial condition Rq 
while Ri is computed as Ri ~ QPRq^Q^^. 

A sufficient condition that points generated according to (J14ap form a peri- 
odic trajectory (a finite set of points) is that for some integer K the following 
conditions hold: 

KP = 2TO7r, 2Ka = 2mr, (16) 

for some integer m, n; in this case the parameter /3/ {2a) = m/n is rational. 

Let us note that due to properties of the Pauli matrices the vector x (7) 
defined by 

a ■ X (7) — S (T ■ X S^^ — exp {i-^cr • s) ct • x exp {—i^u ■ s) , (17) 

where x = [x^, x^, x^\ , s = [s^, s^, s^\ are unit vectors and 7 is a corresponding 
angle of rotation, is easily computed as 

X (7) — cos (7) X + sin (7) X X s + (1 — cos (7)) (s • x) s. (18) 



Using this result it is possible to write the solution (|14ap in a closed form. 
Indeed, equation (J14ap can be explicitly written as 

c-r2if — exp [iKaa ■ q) exp iiK^a ■'p\ ct-Tq exp ( —iK^a ■ p j exp {—iKaa ■ q) . 

(19) 
Applying Eq. (fT8|) twice we find the closed form solution of Eq. (fT3|) : 

r {Ka) = cos {Ka) t [Ka) + sin {Ka) t {Ka) x q + (1 - cos {Ka)) (q • t {Ka)) q, 

(20a) 

t {Ka) — cos {\Ka) Vq + sin {\Ka) Tq x p + (1 — cos {\Ka)) (p • ro) p, (20b) 

where r {Ka) = Y2k, \ = ^I (2a) and i^ = 0, 1, 2, ... . 

The sequence of vectors Tq , r2 , r4 , . . . , generated according to Eq. ((20)) , sam- 
ples a continuous curve C. Indeed, it follows from Eg.f [^D|) that for very small 
a and fixed A vectors r {Ka) and r {{K + 1) a) are very close on a unit sphere. 
Therefore for decreasing a, /? and fixed A = /?/ (2a) the sequence Tq, r2, r4, . . . 
approximates the curve C more and more exactly. It is useful to introduce new 

variable d = Ka which can treated as continuous since a is arbitrary. Using 



this we can rewrite Eq. (|T9)) as: 

(T-r {ff) ~ exp {iQa ■ q) exp {i\Q a ■ p) cr-ro exp {—i\9a ■ p) exp {—iOa ■ q) , (21) 

while Eg. ([20]) leads to explicit formula for the curve C: 

r {9) = cos {9) t {9) + sin {9) t (61) x q + (1 - cos {9)) (q • t {9)) q, (22a) 

t (6*) = cos {\9) ro + sin {\9) ro x p + (1 - cos {\9)) (p • ro) p. (22b) 

3.3 Symmetry and restrictions of dynamics 

Dynamical system (I13|) has continuous symmetry: 

Rn -^ Q'^RnQ^'", Vk e K. (23) 

This symmetry is equivalent to rotation of rjv around q about an arbitrary 
angle Ka. It can be thus expected that dynamics of the quantity rjv ■ q should 
decouple from other degrees of freedom in P^ [US]- Indeed, it follows from 
(PHal that 

r{9)-q = t{9)-q. (24) 

Since |t(6')| = |q| = 1 it follows from the Schwartz inequality that —1 < 
t (^) • q < 1- Now, for given p, q and ro we have 

Ai<t{9)-q<A2, -l<Au A2<1. (25) 

The constants Ai^2 depending on the parameters p, q and the initial condition 
ro can be computed from (j20b[) by elementary means 

Ai.2 = c T ^Jb^ + {a-c)\ (26a) 



where 

a = ro • q, 6 = (ro x p) • q, c = (p • Tq) (p • q) • (26b) 

It thus follows that the motion on the sphere |r (6')| = 1 is bounded by two 
parallels: Ai < r (6*) • q < A2. 

4 Differential equation 



The map p3)) , as follows from the solution shown in Eq. (pOj) , is the stroboscopic 
map with strobe time a of a differential equation which will be deduced from 
the solution (|21l ) where is a continuous variable (it can be derived from the 
form (P^ as well). Differentiating Ea. ipTj) with respect to 6 and using p^ we 
get (see Exercise 41.3 in Ref. [H] for similar computations) 

'^'"'"'^^^ =t[a-u{e), a -rid)], (27) 



d9 
where [A, B]=AB- BA and 



u(e)-q + Ap(0), A=i^, (28a) 

CT • p (0) = exp {i9(T • q) cr • p exp {—iOa ■ q) . (28b) 

It follows that the sequence ro,r2,r4, . . . , generated from Eq. pp)) . samples 
the flow of Eq. (|27|) exactly. As it was remarked in Subsection 3.2 the vectors 
r2K and r2(i^:-|-i) are very close on a unit sphere for very small a and fixed A, 
cf. Eq. ([20| . Therefore for decreasing a, j3 and fixed \ = (3/ (2a) the sequence 
of points on the unit sphere given by rg, r2, r4, . . . can be arbitrarily dense. 

Equations ([27]), (EE]) can be written in exphcit form. Using Eqs.(IT7|), p^ 
we get 

p (61) == cos {26) p + sin {26) p x q + (1 - cos {26)) (q • p) q, (29) 

and hence 

^ = 2r(0)xu(0), (30) 

with u{9) given by Eas. (|28ap . ((29| . Obviously the length of the vector r is a 
conserved quantity, and we shall put \v {6)\ = 1. The angle 6 can be treated as 

do 
increasing with angular velocity w — const, — — lo 121 , and time variable can 

dt 
be introduced to obtain finally 

^-2r(t)xu(i), (31a) 

dt 

u(t) =:a;{[l + A — A cos {2ujt) (q • p)] q + A cos (20-1^) p + A sin {2ujt) p x q} . 

(31b) 

Let us note that Eq.((3T|) is the Bloch equation [Mj. Equation ([31]) has two 
invariants: |r(i)| — const, u (t) • — ^j-^ = which follow from the structure of 

dnnD. 



5 Computational results 

We have performed several computations for the Bloch equation ([31]) and the 
discrete-time dynamical system (1131) . parameterized as described in Section 3.2, 
to show dynamics of the Bloch equation and to demonstrate how the map (|13p 
samples the flow of Ea. ([?T|) . Exact solutions of the map p^ as well as of the 
Bloch equation ([?T|) are given by (^0)) and (HH), respectively. 

Since lo determines time scale only we put uj = 1. In all computations 
described below vectors p, q are orthogonal, p — [1, 0, 0], q = [0, 0, 1] and 
the initial condition is Yq = [0.6, 0, 0.8]. Motion on the sphere is bounded by 
parallels Ai^2 = Tmaxlr (0) • q| given by (pS)) . For the present choice of p, q 
, ro we have Ai^2 = TO. 8. In all figures below the vector q, parallels Ai,2 and 
the equator are plotted, where thin dashed lines indicate points invisible for the 
observer. 

The solution (PO)) . of discrete-time dynamical system (IT^ with Q, P given by 
(dSl) has been plotted in Fig. 1 for a = 0.01, /3 = 0.06 (A = /3/ (2a) = 3) where 
angles a, /? have been given in degrees. The value of a is so small that points 
ro, r2, r4, ... lie so close one to another that a seemingly continuous curve, 
sampling the flow of the Bloch equation ((3T|) very densely, has been obtained. 
The dynamics has been also generated directly from (fTS]) with R\ — QPRq^Q~^ 
(the angle Xo 7^ and arbitrary otherwise) to the same effect, the unit vectors 
Fat have been renormalized after each iteration to avoid numerical instabilities. 
The whole trajectory has three-fold symmetry with respect to the q axis. 




Figure 1: The Bloch equation ((3T|) and discrete-time dynamical system ([13 
a = 0.01, /3 = 0.06, A = 3. 



In Fig. 2 dynamics of vectors r^ obtained from ((20|) has been plotted for 
a — 2, /3^8{X^/3/ (2a) = 2). We thus obtain forty five points marked with 
dots. The solution ^ for a = 0.01, (3 = 0.04 (A = /3/ (2a) = 2) sampling the 



Bloch equation densely has been also plotted. The closed curve has two-fold 
symmetry with respect to the q axis. 




Figure 2: The Bloch equation ([3T|) and discrete-time dynamical system (|T3 
a = 2, /3 = 8 (dots) and a = 0.01, /? = 0.04 (sohd and dashed lines), A = 2. 




Figure 3: The Bloch equation ([3T|) and discrete-time dynamical system ([T3|) . 
a = 0.01, /? = 0.0205, A= 1.025. 

In Fig. 3 initial stage of dynamics of vectors rjv ((20l ) has been plotted for 
a = 0.01, (3 = 0.0205 (A = (3/ (2a) = 1.025), dot on the upper parallel marks 
the initial vector rp. Had the value of A be equal to one exactly the trajectory 



would consist of one closed loop. Since A — 1.025 = 41/40 the whole closed 
curve consists of forty one loops (in Fig. 3 seven such loops have been shown) 
and samples densely the Bloch equation. In the case of close but irrational value 
of A, e.g. A — \/2 — 0.389 = 1.0252 . . . , the trajectory is ergodic on the whole 
spherical sector bounded by two parallels A\_2 ~ =f0.8. 

6 Summary and discussion 

We have introduced in Subsection 3 a class of discrete-time invertible maps 
([7]) on an arbitrary group Q and the exact solution (fTO|) of this map has been 
found. Maps of form ([7]), parameterized on a Lie group, generate points in 
the parameter space which sample a trajectory in this space. This curve can 
be generated forward as well as backward from a given initial condition. This 
suggests that the group action ([7]) may correspond to a flow of a differential 
equation. The results described above generalize significantly our earlier findings 

m- 

In Section 3.2 a simplified dynamics (fT3|) has been parameterized on the 
unit sphere, i.e. for Q — SU (2), and it has been shown directly by constructing 
the solution ()14a|) . ([20|l and ((22|) that the maps samples the fiow of the Bloch 
equation ([3T|) exactly and arbitrarily densely. It can be thus stated that the 
Bloch equation (PT|) is a unique suspension of the map P^ . 

In the special case p • q = we recover the solution obtained by H.K. Kim 
and S.P. Kim [22], see also Kobayashi papers where several methods to solve 
the Bloch equation were described earlier [23l [24l [25] . It seems that the present 
formulation leads to some progress in understanding the Bloch equation. The 
role of the parameter A has been elucidated in Subsection 3.2. More exactly, 
it follows from the condition P^ that for fixed rational value of parameter 
A = /3/ (2a) the fiow of Ea. ([3T]) generates closed trajectories consisting of a finite 
set of points. For a, /? ^ and fixed value of A points obtained according to 
([20]) sample trajectories of ([31]) exactly and arbitrarily densely. Moreover, it has 
been shown in Section 3.3 that dynamical system (I13p has rotational symmetry 
with respect to the q axis from which dynamical restrictions for the quantity 
r (0) • q have been derived, see ([26]) . It follows from ([22]) that for A = m/n, with 
m, n relatively prime, the curve r {9) has m - fold symmetry with respect to the 
q axis. All these results have been visualized in Section 5 where computational 
results have been presented. The figures have been produced by code written 
in the MetaPost picture-drawing language [25]. 

The results described in the present paper can be generalized in several di- 
rections. First of all it should be determined when the general group dynamical 
system ([7]) samples a continuous curve. Whenever this is the case it should be 
possible to construct from a solution (llOp . computed for some Qn's, an ODE 
for which the group action is a stroboscopic map. Of course, uniqueness of such 
construction should be investigated. A very simple choice of Qn's is Q2K = S, 
Q2K+1 — T, K — 0,1,... . We found in our early computations that for 
ST 7^ TS dynamics of Eq.([7]) was very complicated [27], yet there is a closed 



form solution plT]) . Finally, group actions on other groups such as SL2 (R) can 
be considered. 
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